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ABSTRACT 

We study the wave propagation modes in the relativistic streaming pair plasma of the 
magnetospheres of pulsars and magnetars, focusing on the effect of vacuum polariza- 
tion. We show that the combined plasma and vacuum polarization effects give rise to a 
vacuum resonance, where "avoided mode crossing" occurs between the extraordinary 
, mode and the (superluminous) ordinary mode. When a photon propagates from the 

vacuum-polarization-dominated region at small radii to the plasma-dominated region 
at large radii, its polarization state may undergo significant change across the vacuum 
resonance. We map out the parameter regimes (e.g., field strength, plasma density 
and Lorentz factor) under which the vacuum resonance occurs and examine how wave 
i-S^ ' propagation is affected by the resonance. Some possible applications of our results are 

discussed, including high-frequency radio emission from pulsars and possibly magne- 
Q ' tars, and optical/IR emission from neutron star surfaces and inner magnetospheres. 
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1 INTRODUCTION 



?H ' The magnetospheres of pulsars and magnetars consist of relativistic electron-positron pair plasmas, plus possibly a small 
5^ I amount of ions. These plasmas can affect the radiation produced in the inner region of the magnetosphere or the stellar surface. 
Understanding the property of wave propagation in pulsar/magnetar magnetospheres is important for the interpretation of 
various observations of these objects. 

Radio emission from pulsars (at least "normal"-non millisecond - pulsars) likely originates from close to the stellar surface, 
within a few percent of light cylinder radius (e.g., Blaskiewicz et al. 1991, Kramer et al. 1997). A number of studies have 
been devoted to the propagation effect of radio waves in pulsar magnetospheres (e.g., Cheng & Ruderman 1979; Barnard & 
Arons 1986; Barnard 1986; Lyubarskii & Petrova 1998; Melrose & Luo 2004; Petrova 2006). Some of the observed polarization 
properties of pulsar emission, such as orthogonal modes (in which the polarization position angle exhibits a sudden ~ 90° 
jumps; e.g. Stinebring et al. 1984a, 1984b) and circular polarization (e.g., Radhakrishnan & Rankin 1990; Han et al. 1998; 
You & Han 2006) may be explained by the propagation effect. In addition, optical/IR radiation may be produced in the inner 
magnetosphere or surface of magnetized neutron stars. For example, while for most radio pulsars the optical and near IR flux 
is thought to be dominated by magnetospheric emission, several middle-aged pulsars (PSR B0656-f 14, PSR B0950-I-08 and 
Geminga) also exhibit a surface optical component (e.g., Mignani, de Luca & Caraveo 2004; Kargaltsev et al. 2005, Mignani et 
al. 2006). The optical emission detected in several thermally emitting, isolated neutron stars mostly likely has a surface origin 
(e.g., Kaplan et al. 2003; Haberl et al. 2004; Haberl 2005; van KerKwijk & Kaplan 2006). Finally, the optical/IR emission 
detected from a number of magnetars may originate from a hot corona near the stellar surface (Beloborodov & Thompson 
2006). 

Wave modes in pulsar magnetospheres have been studied in a number of papers under different assumptions about 
the plasma composition and the velocity distribution of electron-position pairs (e.g., Melrose & Stoneham 1977; Arons & 
Barnard 1986; Lyutikov 1998; Melrose et al. 1999; Asseo & Riazuelo 2000). In this paper, we reinvestigate the property of 
wave propagation in the magnetospheres of pulsars and magnetars, focusing on the competition between the plasma effect 
and the effect of vacuum polarization. It is well known that in the strong magnetic field typically found on a neutron star. 
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the electromagnetic dispersion relation is dominated by vacuum polarization (a prediction of quantum electrodynamics; e.g., 
Heisenberg & Euler 1936; Adler 1971; see Schubert 2000 for extensive bibliography) at high photon frequencies (e.g. X-rays) 
and by the plasma effect at sufficiently low frequencies (e.g., radio waves). But where is the "boundary" at which the two 
effects are "equal" and what are the mode proportion in the "boundary" regime? These are the questions we are trying to 
address in this paper. We show that the combined plasma and vacuum polarization effects give rise to a vacuum resonance: 
For a given plasma parameters and external magnetic field, there exists a special photon frequency at which the plasma effect 
and vacuum polarization effects "cancel" each other. A more physical way to describe the resonance is as follows: Consider a 
photon propagating in the inhomogeneous pulsar/magnetar magnetosphere with varying plasma density (and/or distribution 
function) and magnetic field. For certain parameter regimes of the photon frequency, plasma density and magnetic field 
strength — to be determined in the following sections, the photon may traverse from the vacuum-polarization-dominated 
region to the plasma-dominated region or vice versa. This transition point (location) is the vacuum resonance. When the 
photon crosses this resonance, its polarization state may undergo significant change. The goal of our paper is to map out the 
parameter regimes under which the vacuum resonance may occur and to elucidate how wave propagation may be affected by 
the resonance. 

Vacuum resonance in cold, non-streaming plasmas have been studied before (e.g., Gnedin et al. 1978; Meszaros & Voutura 
1979; Lai & Ho 2002, 2003a). In the atmospheres of highly magnetized neutron stars, the resonance can significantly affect 
the surface emission spectrum and polarization (Ho & Lai 2003; Lai & Ho 2003a, b; van Adelsberg & Lai 2006). We note that 
while some previous papers on wave modes in pulsar magnetospheres (e.g. Arons & Barnard 1986) did include the vacuum 
polarization contributions to the dielectric tensor, the vacuum resonance phenomenon was neglected because it is unimportant 
at the low frequencies and magnetic fields they considered. 

Our paper is organized as follows. In §2, we give the expression for the dielectric tensor of a relativistic pair plasma 
charactering the magnetosphere of pulsars/magnetars, including the contribution due to vacuum polarization. In §3, we 
derive the general expression for wave modes in the combined "plasma -|- vacuum" medium, and show that the vacuum 
resonance arises for a wide range of magnetosphere parameters. In §4 we study the evolution of wave mode across the vacuum 
resonance. In most of this paper, we consider cold, streaming plasma with a single Lorentz factor 7. We examine the effect 
of a more general 7 distribution in §5 and the case of opposite plasma streams in §6. In §7 we discuss possible applications of 
our results. 



2 DIELECTRIC TENSOR FOR AN STREAMING ELECTRON-POSITRON PLASMA 

Wo consider an electron-positron plasma in the magnetosphere of a neutron star (NS). Let A'^_, A'^+ be the number densities 
of electrons and positrons, respectively, N = N- + N+ the total density, and / = A''+ /N the positron fraction. The corotation 
region of the magnetosphere is usually assumed to have Goldreich-Julian charge density 

^= = -2^"--^' (2-^) 
where 17 is the angular velocity of the star. This is not necessarily the case in the open-field line region. In this paper we shall 
use the Goldreich-Julian number density as a fiducial value: 
O R 

Ng3 = ^ 7.0 X 10'°Si2Pr' cm-^ (2.2) 

2-Kec 

where B12 = 5/ (10^^ G), Pi is the spin period in units of 1 s. The actual particle density N is larger than A^gj by a factor of 
T], i.e. N = tjATgj, with > 1. If the charge density is equal to the Goldreich-Julian value, then r]{l — 2/) = 1, but we will not 
entirely restrict ourselves to this constraint. 

Although the plasma is expected to be relativistic, it is useful to define the (nonrelativistic) cyclotron frequencies and 
plasma frequencies of electron and position: 

eB 

UJc = U)c+ = U!c- ~ , (2.3) 

rrieC 

= (1 - f)ujp, (2.4) 



4nN+e^ , 2 



= H, (2.5) 



2 477 jVe^ 

= (2.6) 

The characteristic values are 

z/c = ^ = 2.795 X 10'' B12 GHz (2.7) 
27r 
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;/p = 5^ = 8.960 X 10^ Af^/^ Hz = 2.370 r]^^^Bl^''P~^^'^ GHz, (2.8) 
27r 

Because of the very short cyclotron/synchrotron decay time of elections and positions (« 3 x 10~^®-B]^^7 sec) , all 
the particles in the magnetosphere quickly lose their transverse momenta and stay in the lowest Landau level. Thus the 
magnetosphere pair plasma can be considered as one-dimensional, with the particles streaming along the field line. The Lorentz 
factor 7 of the streaming motion is uncertain. In the polar-cap region of a pulsar, primary particles may bo accelerated to 
very high energy (7 ^ 10® — 10^) by a field-aligned electric field. The bulk of the plasma produced in an electromagnetic 
cascade may have lower energies, 7 ~ 10^ — 10'', with multiplicity factor 77 ~ 10^ — 10® (e.g., Dangherty & Harding 1982; 
Hibschman & Arons 2001). Kunzl et al. (1998) argued that too high a density of secondary particles in the magnetosphere 
is in contradiction to the observed low-frequency emission from radio pulsars, implying 77 ^ 100. The physical parameters for 
the plasma in the closed-field-line region of a pulsar are also not well constrained. It was suggested that a pair plasma density 
larger than Nqj may be present, maintained by conversion of 7-rays from the pulsar's polar-cap and/or out-gap accelerators 
(see Wang et al. 1998; Ruderman 2003). 

For magnetars, recent theoretical work suggests that a corona consisting mainly of relativistic pairs with 7 ~ 10^ (and 
a wide spread in 7) may be generated by crustal magnetic field twisting/shearing due to starquakes (Thompson et al. 2002; 
Beloborodov & Thompson 2006). The plasma density is of order A?" |V x B|/(47re) B/{4Trer) (for a twist angle of order 
unity), implying r] = N/Nqj ~ c/{2Q,r) ~ 2 x 10*(Ji/r) (where R is the stellar radius). There are roughly equal amount of 
electrons and positrons, / ~ 1/2, with the electrons and positrons streams in opposite directions. 



2.1 Cold, Streaming Pair Plasma 



We first consider a cold electron-positron plasma with all the electron streaming with velocity V_,o and positron all with 
V +fi which is also along the magnetic filed Bo- The dielectric tensor for such a plasma was derived by Melrose & Stoneham 
(1977) based on Lorentz transformation (see also Melrose 1973). Here we outline a derivation based on classical magneto-ionic 
theory. 

The equation of motion of a given charge species (mass rus = me, charge qg = ±e, s = ±) reads 



^ {ismsVs) = qsE + ^Vs X B, 



(2.9) 



where l3s = VJc, 7^ = (1- K^/c^)"^/^, = Vgfi + SVg, E = SE, B = Bq + SB. Here SVs, SE, SB are associated with the 

disturbance in the plasma, and have the form e*^*" ' With dSVs/dt — —iujdVs+ic {k ■ Vs.o) SVs and SB = (ck/u)) x SE, 
wo can solve for SV s in terms of the components of SE. The current density associated with the disturbance is 

SJ = "^{N.qgSVs + dNgqgVg^o), (2.10) 

s 

where the sum runs over each charged particle species s (electron e and position p). The density perturbation SNs is determined 
by the continuity equation, dNs/dt — —V ■ (NsVs.o), which gives SNs ~ Ns{k ■ SVs)/{i^ — k ■ ^,5,0). The conductivity tensor 
is defined hy SJ = a ■ SE, and the dielectric tensor is given by e = l + i{4nc/uj)a, where I is the unit tensor. In the coordinate 
system x'y'z' with Bo/\Bo\ along z', and fe in the x'-z' plane, such that k x Bq = — sin^sy' {6b is the angle between fe and 
Bo), we find 



S iD 
-iD S 
A iC 



A 

-iC 
P 



(2.11) 



with 

S = 1 + 



D 



s 
s 

P = 1 + I](^'7+C'/mi), 



(2.12) 
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where 

fs.ll 

/a, 12 
fs,ri 



1 — Usjs ^ (1 - nl3s COS Ob) ^ ' 
sign (qs) ul^^Vs'yr^ (1 - nps cos6>b)~^ 
1 - Us-yT^ (1 - nf3s cos 6s)~ 

Vs 



7J (1 - n/3s cos 6b) 



2 ' 



C. = , . (2.13) 

1 — n/3s cos Ob 

(the subscript "0" has been suppressed), n = ck/uj is the refractive index and the relevant dimensionless quantities are 
2 

u-=u+=u='^, (2.14) 



2 2 2 

= «+ = ^, « = ^ = 2.15 

oj^ oj'^ a;'' 



a.i.i Pair Plasma With the Same Velocity (P- = (3+ = (3) 

Consider the case where all the electrons and positrons have the same velocity {V - = V+, f}- = P+ = (3, 7_ = 7+ = 7). 
In this paper, we will focus on the regime where ujc 3> "fto (1 — n/3cosSs), or U7~^ (1 — nf3 cos6b)~^ 3> 1, i.e. the photon 
frequency shifted to the plasma rest frame is much lower than the cyclotron frequency. In this regime, the components of 
dielectric tensor can be simplified to 

S = I + /11, 

D = fl2, 

A = C/11, 
C = C/12, 

P = l + ,A+CVii, (2.16) 
with 

fii = ~ vu~^j (1 - nf3 cos 9b)^ , 

S 

fi2 = ^/a,i2i:-(l-2/)i;w"'''^(l-n/3cos6lB), 

S 

fn = fs,r, ^ -f7~^ (1 - nf3 cos6b)~^ , 

S 

C = n/3sin6ls (1 -n/3cos6ls)"^ (2.17) 
2.1.2 Pair Plasma with opposite velocity (the case of — —l3+ = P) 

Suppose the plasma is composed of two opposite streams: one is the electron stream with /3_ = /? and the other is the positron 
stream with /3+ = — /3. For ojc ^ 70^(1 ± nP cosOb), Eqs. (|2.13|) simplify to: 

~ fvu~^^ {1 + nPcosOBf , 

f-^ii ~ (1 - /)i;it"^7(l - n/?cos6's)^ , 

/+_i2 ~ fvu^^^^ {1 + uPcosOb) , 

f-,12 ^ ^{l- f)vu-^^^ (l-nPcoseB), 

.f+,v - -fvy^'^ {I + nP cose b)"'^ , 

f-.v - ^{1^ f)v-f-H^~riP cose Br\ 

^± = =Fn/3sin6ls (1 ±n/3cos0s)"V (2.18) 

2.2 Pair Plasma with a Distribution of 7 

The pair plasma in pulsar magnetosphere many have some spread in the Lorentz factors, although the precise distribution is 
unknown. For a general distribution function fsija), normalized J fsi'fs) d'ye — 1, we can average our result in §2.1 to obtain 
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the dielectric tensor components: 

s 

s 
s 
s 



(2.19) 



These expressions agree with the result derived using standard kinetic theory [e.g. Krall & Trivelpiece 1986, equation (8.10.11); 
Arons & Barnard 1986; Melrose & Stoneham 1977, Lyutikov 1998]. 

Given the uncertainty in the 7-distribution, we shall consider the simplest flat distribution, for both electrons and 
positrons: 

;(^)^/ V(2A7) 7c-A7<7<7c + A7 

\ otherwise ^ ' 

where we shall assume 7min = 7c — A7 ^ 1 (so that |/3| ~ 1 for all particles) and (7c + A7)lj(1 ± nP cos 9b) i^c (so that 
the Doppler-shifted photon frequency is always below cyclotron resonance). Then the components of the dielectric tensor can 
be evaluated analytically. 

For the case of = /3+ = /3, substituting Eqs. (t2T7)l and (p:20]l in Eq. (f2T9)l . we obtain 

S = 1 + Fii, 
D = F12, 
A = CFii, 
C = CFi2, 

P = 1 + F, + C'-Fii, (2.21) 
with 

Fii ~ vu~'^'Yc {1 ~nP cos Ob f , (2.22) 

F12 ^ -{l-2f)vu-^^^{l~nPcoseB), (2.23) 

Fr, = -7;7-^ (I-A7V7C)"'' (1-^/3 cos efl)-^ (2.24) 

Note that Fu and F12 are unchanged compared to the case of delta- function distribution (Eg. I2.17)l . while is changed by 
a factor (l - A7^/7^)~^ 

For the case of jS- = -/3+ = /3, substituting Eqs. ([2TT3)) and in Eq. (1271911 . we obtain 

S = l + ^F,,ii, 

s 

^ = ^CsFs,ll, 

s 

s 

P = 1 + X(F,,„ + C?F,,ii), (2.25) 

s 

with 

^ fvu-\c il + nl3 cos Ob f , (2.26) 
^ (l-/)TO"Sc(l-n/3cosfe)^ (2.27) 
F+,12 =i /uit~'/^ (l + n/3cos6lfl), (2.28) 
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F+,n = -/«7c"'(l-A7V7?)"'(l + n/3cos0s)-\ 
F-,r, = -(l-/)^7-^ (l- A7V7c)~'(l-n/3cos0i3)- 



(2.29) 
(2.30) 
(2.31) 



2.3 Correction due to Vacuum Polarization 

Vacuum polarization contributes a correction to the dielectric tensor: 
Ae'^> = (a - l)l + qi3B, 



(2.32) 



where I is the unit tensor and B = B/B is the unit vector along B (here we use B, B to denote Bq, Bq for simple notations). 
The magnetic permeability tensor fi also deviates from unity because of vacuum polarization, with the inverse permeability 
given by 

fj.^ = al + mBB. (2.33) 

In the low frequency limit Tiuj <^ rricC? , general expressions for the vacuum polarization coefficients a, q, and m are given in 
Adler (1971) and Heyl & Hernquist (1997). For B <^ Bq = mlc^/{eh) = 4.414 x lO" G, they are given by 



a = 1 - 2Sv, 
where 

X ap f B 
^^^45^1b^ 



7(5v 



-45v 



= 2.650 X 10~*B?2 



(2.34) 



(2.35) 



and dp = /Tic — 1/137 is the fine structure constant. For B ^ Bq, simple expressions for a, q, m are given in Ho & 
Lai (2003) (see also Potekhin et al. 2004 for general fitting formulae). 

When |Ae,'j>| < 1 or B/Bq <C Stt/qf (B < 5 x 10^® G), the plasma and vacuum contributions to the dielectric tensor 
can be added linearly, i.e., e = e'''^ + Ae^^'. In the frame with B along z , 



S' iD A 
-iD S' -iC 
A iC P' 



(2.36) 



with S' = S + a, P'^P + a + q and a = a - 1, 



3 WAVE MODES IN A COLD STREAMING PLASMA 

Here we consider the case of a pair plasma all streaming with the same velocity /3 along the field line. The effect of finite 
spread in 7 will be studied in 351 and the case of opposite streams will be considered in §6. 



3.1 Equations for the Wave Modes 

Using the electric displacement D = e ■ E and equation (|2.33|) in the Maxwell equations, we obtain the equation for plane 
waves with E oc e^f*' '--'^') (henceforth we use E to denote SE, and use B to denote Bo) 

1 



r 1 2 

la 



-(k X B)i(k X B)j 



E, = 0, 



(3.37) 



where n = ck/to is the refractive index and k — k/k. In the coordinate system xyz with k along the z-axis and B in the x-z 
plane, such that k x B — — sin^sy, the components of dielectric tensor are given by [compared to eq. (|2.36p ] 

S' cos^ Ob - 2 A sin 61b cos Og + P' sin^ Ob , 
S', 

S' sin^ 6b + 2 A sin 6b cos dB + P' cos^ 6b, 
~<^yx = i {Dcos6b — C sin 6^3) , 

^ Acos2eB + {S' - P')smeBCOs6B, 
—e^y = ~i{D sin6B + C cosdB) , 
The z-component of equation (|3.37l) gives 
Ez — —£3/ {tzxEx + ezyEy) . 



^yy 

^xy 
^xz 
^vz 



(3.38) 
(3.39) 
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Reinserting this back into equation (|3.37|) yields 

2 



Vxx — n 



Tjxy 



r]yx Tjyy TTl 



Vyy 



Vyx 



Ex 

Ey 
,2 , 



= 0, 



where r = l + f = l + (m/a) sin 6b and 



i^zz^xx ^xz^zx ) iSP A 

a^zz Ci^zz ^ 



(^^zz^yy ^yz^z 



— [(S'^ - - S'P' + C^) sin^ Ob + 2(AS' - CD) sin^s cos6ifl + S' P' - C^l 



■ {^zz^y 



V^y y^zz^yx '^yz^zx ! 



— - {P'OcosOb - S'CsmeB +A(DsmeB - C cos 6is)] 



(3.40) 



(3.41) 



(3.42) 



(3.43) 



3.2 B = oo limit without QED effect 



The above expressions are valid for the general dielectric tensor f Eg. 12.36]) . We now consider the B = oo limit introduced by 
e.g., Tsytovitch & Kaplan (1972) and Arons & Barnard (1986). In this regime, the magnetic field is sufficiently large so that 
the cyclotron frequency are large compared to the Lorentz-shifted wave frequency. At the same time, B is not really infinity 
so that the wave propagation is dominated by plasma effect and we neglect the QED correction in the dielectric tensor. The 
approximate elements of the dielectric tensor are 



/ii ^ vu ^7 ~ 0, fi2 ~ —(1 — 2f)vu 
fr, ~ ~vy^''* (1 - uPcosObY^ , 



-1/2 



0, 



S^l, P~l + /„, -D~0, A~0, C~0, 

Ezz =i 1 + /r, COS^ 6b ■ 

Equations (|3.41|l - (|3.43|l then reduce to 



fjxx 



1 + /., 



1 + fr, COs2 9b ' 

Vyy — 1 1 

Vyx — Vxy — 0- 

Solving equation (|3.40|) . we find (see Arons & Barnard 1986) 



^ — Vxx 



1 + /. 



1 + fn cos^ 



(3.44) 
(3.45) 
(3.46) 
(3.47) 



(3.48) 



(3.49) 



(t^' - c^fcf ) 



c^fci = 



73a;2(l -/3cfc||/a;)2 
(where A;|| — fccosSs, ki_ = fcsin6'_B) for the ordinary mode (0-mode) and 



2 1 

n = 1 



for the extraordinary mode (X-mode). The polarization of the O-mode is given by 







Vyx 


~0, 






Ex 




Vyy ~ VxX 






E± 




fr, sin 9b cos 9b 






Ex 




1 + cos 


^Bb 




tan 



and that for the X mode is 
0. 



Ex 




El 


E~y 


- 0, 


Ey 



(3.50) 

(3.51) 

(3.52) 
(3.53) 

(3.54) 
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Figure 1. The refractive indices n = ck /w of wave modes as a function of LUp/w for a cold plasma with 7 = 1 and the wave propagation 
angle 8b = 30°. The solid line shows the superluminous O-mode, the dashed line the X-mode, and the dot-dashed line the subluminous 
O-mode. The insert shows the blowup of the "avoided crossing" region between the X-mode and the superluminous O-mode due to 
vacuum resonance. The mixed modes are labeled mode and "— " mode. For the insert, the x-axis gives lO'^Wp/o;, and y-axis gives 
10* (n - 1), and the other parameters are B = 10^^ G, N = Nqj with P = 1 s. 




Figure 2. Same as Fig.l except for 7 = 1.1, and the x-axis gives {ujp/uj)^ "^1"^. Note that for (ujpl'S)'^ '^/^ <g; 1, the subluminous 
O-mode has n = ~{u)p /u))^~^/^ + l/{f3cos9g), and it is shown as a dotted line in the figure. 



Thus, the X-mode is a transverse wave with the electric field vector in the k x B direction, while the O-mode is polarized in 
the plane spanned by k and B. 

Figures [1] - [3] depict the refractive indices n = ck/ui of different modes as a function of {ujp/uj)'y^^^'^ for 7 = 1, 1.1 and 
10"^, respectively, all with 6b = 30°. The O-modes have two branches: the superluminous branch {w > ck or n < 1) and the 
subluminous branch {w < ck or n > 1), the latter corresponds to plasma oscillations. At low densities, {ujp/u!)'y~'^^^ <^ 1, the 
superluminous O-mode becomes transverse vacuum electromagnetic wave which can escape from the magnetosphere. 

In the very low density region, (a;p/tj)7"^/^ < 1, the QED effect may not be neglected compared to the plasma effect. 
The "competition" between the vacuum polarization effect and the plasma effect gives rise to a vacuum resonance, at which 
the superluminous O-mode and the X-mode may be coupled with each other. In the remainder of this paper, we will focus 
on this vacuum resonance phenomenon. 
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Figure 3. Same as Figure [2l except for 7 = 10^. In the insert, the x-axis gives 10^(tL'p/a;)7 "^/^ and the y-axis 10*(n — 1). 



3.3 Wave Modes Including QED Effect 

We now consider iiow vacuum polarization affects the X-mode and (superluminous) 0-mode. Because of tlie possibility of 
"mode crossing" , we label the modes as "+ mode" and "— mode" . We write the mode polarization vector as E± = E±t+E±zZ, 
with the transverse part 

1 



E 



' (l+i^2)l/2 

where iK± = Ex/Ey. Eliminating from equation p.40p . we obtain 
where the polarization parameter /3p is given by 



(3.55) 



(3.56) 



/3p 



(5"2 _ 1)2 _ g,p, ^ (^2^ gjj^2 ^ ^^2 _ ^2j _ jj^-^ sin 261s + {S'P' - A^) (1 



2[P'DcoseB - S'CsmOB + AiDsinOB - Ccosfe)] 
The index of refraction n± of the two modes can be obtained from equation (|3.40p . giving 

r r 

We now consider the case of cold streaming plasma with /3_ 
following conditions: 

«7~^ (1 - /3cos6ifl)~^ = 7.812 x 10^^ B?2;^r^7:^^ (1 ~ Pcos9b)~'^ > 1, 
vy-^ = 5.617 X 10-='j7Bi2PrVr^73"^ < 1, 



(3.57) 



(3.58) 



P+ = l3. We focus on parameter regimes satisfying the 



(3.59) 
(3.60) 



where u is the wave (photon) frequency, i^i = !//(lGHz) and 73 — 7/10"^. The first condition implies that the Doppler- 
shifted frequency is lower than Wc, so that we can use equation (|2.17|) for /n, /12 and /,,. The second condition implies 
I/12P <S <S 1, and l/^l <C 1. Under these conditions, equation (|3.57p reduces to 



/3p 



{fv + 1 + ^) sin^ 6b — /ii [1 — (cos 6^3 — (siiiOb)^] 



2/12 (cos f 



sin 9e 



We shall see that conditions p.59p and (|3.60|l also imply that the index of refraction is close to unity. 



(3.61) 



11 < 1. Thus 



(cos t 



C sin 9e 



{cosOb — nf3) /{I — nf3cos0B) ~ 1 — 7 sin^s (1 — P cosOb) , and we can easily check that the 



second term in the numerator of equation (|3.6ip is much smaller than Equation (|3.6ip therefore simplifies to 



fn + q + m 



2/1; 



(cos Ob — C sin 9b) 
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fv sin^ 9b ( ^ ^ q + m ^ 



2/12 (cos 61s - C sin 6*5) \^ j 
= PoPv, (3.62) 
where /3o is the polarization parameter in the absence of vacuum polarization: 



fv sin^ ' 



^ _j_n 

~ 2/12 (cos6lfl - (sm6B) 

- -Ul-2fr^u^^\-'il-PcoseB)- \ „ (3.63) 

2 (cos fcifl — (, sm y_B) 

and /3v is the correction factor due to vacuum polarization: 



fl ~ 1 . 1 + "^ 



1 



^ 1 1+^ _. (3.64) 

(1 — /3cOS 6b) 

Equation p.58|) for the index of refraction n± of the two modes also simplifies to 

n± ~ 1 — msin^ 6b + fii + /12 (cos^b — ("sin^s) K±. (3.65) 
For positive /3p with |/3p| ^ 1, we obtain simple expressions of refractive indices 

= 1 — msin^ 6b + fii + {f-n + Q + in-) sin^ 6b, = 1 — msin^ 6b + /ii, (3.66) 
while for negative /3p with |/3p| S> 1, we have 

n\ = 1 — msin^ 6b + ./ii, n?. = 1 — msin^ 6*5 + /n + (/j? + g + m) sin^ 6b- (3.67) 
For |/3p| = 0, we have K± = ±1 and 

n\ = 1 — msiv? 9b + i /12 (cos^b — Csin^s) . (3.68) 

From equations (|3.66tf3l68|) . we see that n is indeed close to unity when equations (|3.59|) and (|3.60|l are satisfied. 

3.4 Vacuum Resonance 

For |/3p| 2> 1, the two modes are (almost) linearly polarized: the mode with \K\ ~ 2|/3p| 2> 1 is polarized in the k-B plane, 
and is usually called ordinary mode (0-mode); the mode with \K\ ~ l/(2l/3p|) ^ 1 is polarized perpendicular to the k-B 
plane, and is called extraordinary mode (X-mode). From equation (|3.62|) . we see that for a general 6b which is not too close 
to 0° or 180°, and for almost all values of B, N, 7's, the inequality [/3p| 2> 1 is satisfied either when the condition 



fv -in 

—!-sm 63 

/12 



(1 - 2fy^u^^^'y-^{l ~ cos^s)"^ sin^ 6b 



= 2.795Bi2!^rS:r^(l-cos6lB)-^sin^6ls(l-2/)-^ > 1 (3.69) 
is satisfied, or when 



q + m . 2n 

sm Ob 











15^ 





-( (l-2/)-'M'/''t;-'(l-cos6'B)-'sin''ej; 
'TT \Bq J 



= 39.56Bi2i'i (1 - 2/)"^77"^Pi(l -cos6is)"^sin^6is > 1 (3.70) 
is satisfied. The exception occurs when 

fr,+q + m = Q (3.71) 
or 

/3v = 1 "l^^ — = 0. (3.72) 

U7~3 (1 _ I3cos6b) 

This defines the "vacuum resonance". For given u, 7 and B, the resonance occurs at the density 

Nv = 9.905 X 10"B?2i^?73 (1 -/3cos6lB)^F(6)cm"^ (3.73) 



where 

q + rn 
QpV(157r 



11 



o 

T 



oo 
O 



oq 

Csj 



CD 
C\j 

LT) 
CN 




V 73 

Figure 4. The polarization ellipticity K (upper panel) and index of refraction n (lower panel) of the wave mode as a function of the 
plasma density parameter rj = N/Nqj (see equation 12.21 1 near the vacuum resonance for B = 10^'^ G, !^ = 1 GHz, / = and 9g = 45°. 
Three values of 7 are considered: 73 = 7/10^ = 0.5, 1, 1.5. The solid lines are for the "+" mode and dashed lines for the "— " mode. 
Both the resonant density (see Eg. I3.75"] l and width (see Eg. I3.77"] l are proportional to 7"^ for 7 3> 1. 



is equal to unity for b = B/Bq <^ 1 and is at most of order a few for B ^ 10^^ G (see Fig. 1 of Ho & Lai 2003). It's obvious 
that Nv tx 7'^ for 7 S> 1 and Ob not too close to 0° (see Figure |4ll. For 7 = 1, equation (|3.73|) agrees with the result of Lai & 
Ho (2002). The dependence of the resonance density A'^v' on 7 and 9b can be easily understood from the cold (non-streaming) 
plasma limit and Lorentz transform. The wave freqency in the plasma rest frame is cu' — ^{uj — I3k\\ ) ~ 'yuj(l — j3 cos 6b)- Note 
that in this frame, the external magnetic field and the dielectric tensor due to vacuum polarization are unchanged. The rest 
frame plasma density at the vacuum resonance is Ny ex. B'^lo''^ F{h), and the corresponding "lab frame" density is Nv = jNy. 
We can rewrite the resonance density in terms of the Goldreich- Julian density 

Nv 



rjv 



Ng: 



= 14.15PiBi2;/?73 (1 



/3cos Ob) F{h). 



(3.75) 



The physical meaning of the resonance is clear: For given u, 7 and B, the dielectric property of the medium is dominated 
by the plasma effect when 3> A''^, while it is dominated by vacuum polarization when ^ Nv', ai N = Nv, the plasma 
effect and vacuum polarization compensate each other, and the wave modes become exactly circular polarized. 

There are other ways to view the vacuum resonance. For example, at given B, 7 and density A'' (or 77), we can define the 
vacuum resonance frequency: 



0.266 [Pr^B^2^-n'y3^ {I ~ f3 cose By ^ GHz 



(3.76) 



Thus, wave modes with v <^ uv are determined by the plasma effect, while those with v ^ i/v are determined by the vacuum 
polarization effect. 

The characteristic width of the resonance region can be estimated by considering |/3p| = 1 as defining the edge of the 
resonance. Since Pv = 1 — Nv /N, we find that the densities at the edges of the resonance are Nv ± AiV, with 



AN 1 1 3 

— ~ — = 0.7157(1 - 2/) 511^171 (1 ~ 13 cos( 

Nv I Pol 



3 (cos6'_g — C sin 6*3) 



(3.77) 



where cos — C sin 6^5 ~ (cos — /3)/(l — /3 cos 6^). 

Figures |4] and [5] show the mode properties near vacuum resonance for different values of 7 and /. It's obvious that the 
resonance density and width scale with 7^ (for 7 2> 1), and the resonance density doesn't change with /, while the resonance 
region becomes narrow when / is close to 0.5. 
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Figure 5. Same as Fig. |4l except for B = 10^^ G, u = IGHz, 7 = 10'', Ob = 45° and different values of the plasma positron fraction: 
/ = 0, 0.2 and 0.49. 



Note that while vacuum resonance can always be located by + q + m = [Eqs. (13.711) and (I3.72|) ] , significant "avoided 
mode crossing" (as illustrated in Figs. |4] and [5] and the inserts of Figs. [T]-|3]) occurs only if equations (|3.69p and (|3.70p are 
satisfied. If these "linear polarization" conditions are not satisfied, the modes will be approximately circular polarized even 
away from the resonance, and no dramatic change in the mode properties takes place around the vacuum resonance (see 
Fig.©. 



4 MODE EVOLUTION ACROSS THE VACUUM RESONANCE 

Consider a photon (or electromagnetic wave) of a given frequency and polarization state propagating in the NS magneto- 
sphere. The magnetosphere is inhomogeneous because of variations in B, N and possibly 7. How does the polarization of 
the photon evolve, particularly as the photon traverses the vacuum resonance region (e.g. from the plasma-dominated region 
to the vacuum-dominated region)? Clearly, if the variations of the magnetosphere parameters (B, A'^, etc.) are sufficiently 
gentle, the polarization state of the photon will evolve adiabatically, i.e. a photon in a definite wave mode will stay in that 
mode. Then Figure |4] and [5] show that across the vacuum resonance, the photon polarization ellipse will rotate by 90°, with 
the mode helicity unchanged. 

To quantify the mode evolution, it is convenient to introduce the "mixing" angle dm via tanS^ = 1/K+, so that 

tan20™=/3-\ (4.78) 

where we have used |r — 1| ^ 1. The transverse eigenvectors of the modes &rG E jj-jp — cos ^777, , sin ^777, ) and E-T = 
(—1 sin cos 6m). Clearly, at the resonance, 9^ = 45°, the X-mode and O-mode are maximally "mixed". 

A general polarized electromagnetic wave with frequency lu traveling in the z-direction can be written as a superposition 
of the two modes: 

E{z) ^ A+{z)E+{z) + A-{z)E^{z), (4.79) 
Note that both A± and E± depend on z. Substituting equation (|4.79p into the wave equation 

2 

V X (/x"^ ■ V X = ■ B, (4.80) 
we obtain the amplitude evolution equations (see Lai & Ho 2002) 
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Figure 6. Same as Fig. g] except for B = lO" G, u = 100 GHz, 1 - 2/ = 0.1, 7 = 10^, Ob = 45°. Tlie vacuum resonance (defined 
hj fn + q + in = 0) occurs at jj = 1.214 X 10^ (the vertical dotted line). In this case, equations l|3.69|l and l|3.70|l are not satisfied, and 
"avoided mode crossing" does not occur near the vacuum resonance. 



(4.81) 



where ' stands for d/dz, Afc = fc+ — In deriving equation (|4.8ip . we have assumed that E±(z) and j4±(z)exp(— i k±dz^ 
vary on a length scale much larger than the photon wavelength, and we have used k+ ~ fc_ and \k'j-/k±\ <C |fc±|. Clearly, 
when \9'^ \ < |Afc/2|, or 



r = 



26' c 



> 1, (4.82) 



the polarization vector will evolve adiabatically (e.g., a photon in the plus-mode will remain in the plus-mode). Using equa- 
tions (faMIl and (gjH), we find 

eL = -lsin^2g^A f^'^-, ^ . (4.83) 
4 /12 (cos 9b -( sm 9b) Pq 

The difference in refractive indices of the two modes is 

_ ^ /i2(cosgB-CsmeB) 

sm 20rn 

Thus equation (|4.82|l becomes 

r = ^ (1 - 2ff , ^ (1 ~Pcos9Bf (£^^£^£,zi£^)^ » 1, (4.85) 

where H = |/3o//3p| specifies the length scale of variation of Pp along the ray. Equation (|4.85p gives the general condition for 
adiabatic mode evolution along the photon path. 

Clearly, the adiabatic condition (|4.82|) or (|4.85|) is most easily violated at the resonance {9m = 45°). Evaluating equa- 
tion (14.851) at = Nv and using equation (|3.72|l to eliminate cjp, we obtain 

Tv = {f/iye^f, (4.86) 
with the "adiabatic frequency" 
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Figure 7. Same as Fig. |4] except for plasmas witli a flat distribution of Lorentz factor (see equation 12.201 1. The parameters are 
B = 10^2 G, = 1 GHz, / = 0, 7 = 10^, 6*3 = 45°. For each mode, the three lines are for A7/7C = 0, 0.5, 0.8. 



;.ad = 6.081 (1 -2/)-^''S3'' (1-/9 cos es)-^-'''' ( , ) ' R-'^' GRz, (4.87) 

\ cos6'_B — Csinc/s / 

where we have used H ~ \dz/d(3y \ and He — H/{10^ cm). Since /3v = (1 — Nv/N) and /3o changes very slowly near resonance, 
for a constant B, 7, the length H = N/\dN/dz\ (evaluated at A'^ = Ny) becomes the density scale height along the ray. For 
7=1, equations (|4.86|l and (|4.87p agree with the result of Lai & Ho (2002). Using equation (|3.76p to eliminate 7, we can also 
rewrite equation (|4.87|l as 

{l-2f)-''' vt"n-'^' {P^B,,r' F'''H-'/\l~ Pcoseu)-''' ( ^ V^' GHz. (4.88) 

\ cosfc^B — CsmWfl / 



5 THE INFLUENCE OF VELOCITY DISTRIBUTION 

The results of Sections 3 and 4 are for cold streaming plasma with a single Lorentz factor 7. For a general distribution function 
/(7), equations (f537)) - and equations - (I5351l still apply. For the simplest flat distribution of 7 [Eq. (H^^SJ], the 

dielectric tensor functions are given by equation (|2.22p - (|2.24p for 7min ^ 1 and 7max'i'(l — nl3 cos Ob) oJc- If we assume 
1^7";^,-^ <C 1 [see Eq. (|3.60p ]. so that n± ~ 1, we flnd that the vacuum resonance density, width and adiabatic frequency are 
given by; 

riv = — = 14.15PiBi2J^? 7c3(l-A7V7c) "(1-/3 cos 6ifl)"F, (5.89) 

^ = 0.7157 (1 - 2/) B^ij!, (1 - A7V7g)^ d - PcosBsf "^""^"^ , (5.90) 

Jyv ^ ' sm Bb 

v^^ = 6.081 (l-2/)-^/S3-^ (1-/3 cos eB)-V-^/^(' , V^'j/,7^/='GHz. (5.91) 

\ cost^B — Csm&fl y 

Figure [7] shows the polarization ellipticities and the refractive indices of the two wave modes in plasmas with different A7's. 
Thus, a broad distribution of 7's does not qualitatively affect the vacuum resonance behavior. 
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6 VACUUM RESONANCE IN A PLASMA WITH TWO OPPOSITE STREAMS 

We now consider the case where the pair plasma is composed of two opposite streams, for which the dielectric tensor is given 
in Sec. 2. 1.2. The derivation of mode properties for this case are the same as in Section 3. Similar to Eqs. (|3.59p and (|3.60p . 
we assume u')~'^ (1 ± n/3cos S> 1 and v^~^ <^ 1. The polarization parameter l3p fEg. I3.57)l reduces to 

_ (/„,+ + /„,-+g + m)sin2efl + /ii,+ [l-(coses-C+sines)'] - fii.- [l - (cos - C_ sin Ss)'] 
~ 2 [/i2,+ (cos 61b - C+ sin^s) + /12,- (cos^s - C- sin 6(5)] ' 

where we have neglected the terms proportional to /i2,±. We can see that Eq. (|6.92p is rather similar to Eq. (|3.6ip . Since 
cosSs -C±sin6's = = ±1 =F ^7"^ sin^ (1 ± /3 cos 03)"% (6.93) 

1 zt p COS UB 2 

the last two terms in the numerator of Eq. (|6.92p can be neglected compared to the first term. The denominator can be 
simplified to /12.+ - /12,- ^ vu~^^^{l - M/3 cos 6*13), with M = 1 - 2/. Thus, we can rewrite Eq. ((02ll as 

- 2«u-i/2(i_M/3coseB)' ^ ' 



with 



h = ,/.,,+ + /..- - -«7"'— (6.95) 
sm 6b 

and 

Fe = 1 + /3^cos^6Ib +2M/3cose's, (6.96) 

where we have assumed /3 ~ 1, n ~ 1 and sin^s 7^ 0. Similar to Sec. 3. 3, we write /3p — PoPv, where /3o is the polarization 
parameter in the absence of vacuum polarization, and /3v is the correction factor due to vacuum polarization: 

sin^ 9b 



(3o 



2TO-1/2 (1 _ MPcosOb) 

1 1/2 -3 Fe 
— -u 7 



2 ' sin^^s (1 - M/3cos6ls)' 
q + m 

fv 



^ (6.97) 
Vacuum resonance occurs at /3v ~ 0, corresponding to the plasma density (relative to the Goldreich- Julian density) 

= = U.15P,B,,uhf^^Fib). (6.98) 

JVGJ rg 
Similar to Eq. (|3.77|l for the single stream case, the characteristic width of the resonance region is 

AN 3 sin^ 9b {I -MP cos 9b) „„, 
— = 0.7157-Bi2 ^^173 p7 • (6-99) 

Similar to Eq. (|4.85|l . the adiabatic condition is 

cjH 7^^p {1- MP cos 9 
c LUc sin^ {29m) Fg 



r^^ j':l ^i^-MPcosjB)^sinHB^^^ (6.100) 



At the resonance, we have Tv ~ {f/vadf', with 

-1 2/3 



v^i = 6.081 73-^F- 



Fg 



(1 - M/3cos6's)sin'^ 9b 



i?g"'''^GHz, (6.101) 



Figure [8] shows the mode evolution near the vacuum resonance for different values of Afs (or /'s). For the pair plasma 
with M = Q, vacuum resonance still occurs while the resonance region is much wider compared with M > 0. For — 1 < M < 0, 
the modes evolution behaves as the case in Fig. |6] since the linear condition is not satisfied in this parameter regime (given 
in the caption of the figure). 

If there is a velocity spread of the electrons and positrons, for example the flat distribution given by Eq. (|2.20[) . the 
equations above just need to be modified by appropriate factors similar to Eqs. (|5.89p - (|5.91[l : 

m = ^ = 14.15ABi2^?7l(l-A7.V7')'^^n&). (6.102) 

JVG.T rg 
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Figure 8. Same as Fig.|4] except for plasmas with two opposite streams. Note the a;-axis is rjFg / sin^ dg. The parameters are B = 10^^ G, 
u = 1 GHz, 7 = 10^, dg = 45°. For each mode, the three lines are for M = 1 — 2/ = 1, 0, —0.5, respectively. 



= 6.081 7;r'i^(6)-'''^ 



2 sin^ 6Ib(1-M/3cos0b) 



(1 - M/3cosefl)sin^ 6*13 



2/3 



-1/3 , 



(6.103) 
(6.104) 



7 DISCUSSION 

In previous sections, we have studied the property of wave propagation in the magnetospheres of pulsars or magnetars for 
various plasma parameters. We have focused the vacuum resonance phenomenon, arising from the combined effects of plasma 
and vacuum polarization. The possible occurrence of the vacuum resonance and the related wave property depends on the 
plasma parameters, magnetic field and the wave frequency. The key equations are (assuming single-stream plasma): 
(i) The vacuum resonance condition, Eq. p.73p or (|3.75p : 
(u) The adiabatic condition, Eqs. (|4.85I - |488)) : 

In deriving the analytical expressions for the vacuum resonance, we have assumed 

(iii) The Doppler-shifted wave frequency is much less than the electron cyclotron frequency, i.e., uj' = 7cj(l — nPcosOs) ^ t^c 
[Eq. (|339)) ]: 

(iv) The plasma is weakly dispersive, i.e., v^~^ <C 1 [Eq. (|3.60|) ]. 

The vacuum resonance is particularly interesting in the parameter regime such that either 

(v) the waves are linearly polarized away from the vacuum resonance due to the plasma effect [Eq. (|3.69|l ]. or 

(vi) the waves are linearly polarized away from the vacuum resonance due to the vacuum polarization [Eq. (|3.70|l ]. 
Note that at the vacuum resonance, the waves are always circular polarized. 

Figures |9] and [To] summarize these conditions for two different sets of parameters. In both figures, we see that when the 
vacuum resonance induces significant "avoided mode crossing" (cf. Figs. [4] and [5}, i.e. when the resonance lies above the 
"Linear I" or "Linear 11" line, wave evolution across the vacuum resonance is nonadiabatic. In general, this can be understood 
as follows. We define the cross frequency of the "Linear I" line and "Linear 11" line (as well as "Vacuum Resonance" line) as 

Across • 

We find 

i^cross ^ 0.58(1 - 2/)~'''^7?'''^sin4''^Pi"^/S3~^(l - cos6is)"''''^GHz. (7.105) 
Comparing to the adiabatic frequency i^ad, we have 

~ 0.1(1 ~2ff^\^"'p-^/\l ~ cosOb^'^F^/^hI^^ (7.106) 

For typical parameters(e.g., / = - 0.5, 77 = 10^ - 10^, P ~ls, = 45°, He = 1), i-'cross is less than the adiabatic frequency. 
Adiabatic mode evolution across the vacuum resonance with appreciable mode crossing is possible for larger 77 and Ha. 

We now discuss possible implications of our results for various radiation processes in pulsars and magnetars. We assume 
a dipole magnetic field, with 
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Figure 9. Vacuum resonance and wave properties in the magnetic field - frequency domain for the plasma parameters rj = 10, 
1 — 2/ = 0.1, 73 = 1, Og = 45°, Hq = 10 and the NS spin period P = Is. The line labeled "Vacuum Resonance" gives the vacuum 
resonance condition [Eq. II3.75I I]. the line labeled "Adiabatic Frequency" [Eq. I l4.87l l1 means that the adiabatic condition is satisfied to 
the right of the line when the wave evolves across the resonance. The line labeled "a;' < Wc" [Eq. l|3l59)l ] means that the Doppler-shifted 
frequency is much less than the electron cyclotron frequency above the line, and the line labeled "weak dispersion" [Eq. II3.60II ] means 
that the medium is weakly dispersive (i.e., index of refraction close to unity) below the line; Our analytical expression for the wave modes 
and vacuum resonance are valid in this parameter regime (below the "weak dispersion" line and above the "tj' ^ luc" line). Above the 
"Linear I" line [Eq. I I3.69I I] or the "Linear 11" line [Eq. 13.7011 ]. the wave modes are linearly polarized except near the vacuum resonance. 



where B« is the surface magnetic field and Rt the radius of the NS star. Substituting equation (|7.107p into equation (|3.75|) . 
we obtain the location of vacuum resonance (assuming constant rj and 7) 



(TiM-J - (life) m i^r (n) - ^ 



Since the dispersion due to vacuum polarization is of order q + m cc F{b)B^ oc r"", while the plasma effect is measured by 
~ W7~3 oc Ny~^ oc riBj~^ oc rj^'^r'^, if ri^~^ does not vary rapidly, we find that for a given photon frequency u, the wave 
dispersion is dominated by the vacuum effect for r ^ ry and by the plasma effect for r ^ ry- 

First consider the radio emission from the open field line region of a pulsar. The emission angle relative the local magnetic 
field line is 6s ~ 1 /j, so that I — f3 cos 6b — 7~^- This would imply ry /R* <C 1, even for ~ 10^^ G and for high frequencies 
(e.g. V = 20 GHz). Along the ray trajectory, the angle 9b increases. In the small angle approximation {6b 1), we have 

0s^^|y^eo(l-^), (7.109) 

where rem is the radius of the emission point, and do is the polar angle at the stellar surface of the emission field line. Thus 
6b increases from (at r — rem) to (3/4)(rom/ii*)6o. As an example, for rem = 2i?* and 60 ~ \/R*jRhC — 0.0145Pj~^''^, 
equation (|7.1U9|I implies 6b ^ O.Ql^P^^''^ . From Eq. (|7.108p . we find 

This means that for radio emission along open, dipole field lines, plasma effects always dominate the property of wave 
propagation and vacuum resonance will not occur. 

Radio emission may also come from the large-curvature magnetic field structure (e.g., field lines with curvature radius 
~ Rt). In this case, even if 1 at emission, it will become significantly large (~ 45°) after the wave propagates a short 

distance of order R,. Thus, according to Eq. (|7.108p . vacuum resonance can occur for sufficiently high frequencies and strong 
surface magnetic fields. This could be the case with the high-frequency radio emission from the transient AXP XTE J1810-197 
(Camilo et al. 2006). 

Finally, optical/IR radiation emitted from the neutron star surface or near vicinity may experience the vacuum resonance 
while propagating through the magnetosphere. The polarization of such radiation may probe the physical conditions of the 
magnetosphere. 
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